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We analyse the exclusive channel pp — > 7 vr° , assuming handbag dominance. The 
' soft parts are parametrized in terms of CGLN amplitudes for the qq — ► 771"° transition 

■ and form factors for the pp — ► qq ones, the latter represent moments of Generalized 

Distribution Amplitudes. We present a combined fit to Fermilab data from E760 
i-SZh ' taking simultaneously into account information from other exclusive reactions, es- 

pecially from pp — > 77 data. Overall a nicely consistent picture emerges, such that 
one can hope, that our theoretical analysis will be reliable also for the kinematics 
of GSI/FAIR, which hopefully will provide much more precise and complete data. 
Consequently, data from this facility should improve our knowledge both on the 
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proton-antiproton distribution amplitudes and the pion production mechanism. 

Introduction 



The reliable theoretical treatment of hard exclusive processes has been a challenge for 
QCD for many years. With the advent of the generalized parton distribution formalism 
[1] a large class of such processes, all involving some hard scale Q 2 , can now be treated 
on a firm, perturbative QCD basis, absorbing all non-perturbative soft physics in suitable 
generalized parton distributions. This success is made possible by the dominance of the 
handbag diagrams in leading twist, as demonstrated in the factorization proofs. 
There are other hard exclusive processes for which these rigorous proofs do not apply. 
Still, however, there are good arguments for the dominance of the handbag contribution in 
certain regions of phase space also for many of these reactions. Examples of such processes 
are two-photon annihilations into pair of hadrons for large but not asymptotically large 
Mandelstam variables, s, —t, —u, for applications see [2, 3, 4]. 

The FAIR project at GSI with the HESR anti-proton program [5] will offer ideal possi- 
bilities to study exclusive channels in pp annihilation, e.g. pp — > 77, the time reversed of 



the widely studied process 77 — > pp [6, 7]. Another very interesting channel is pp — > 7r°7 
because rates will be much higher and the amplitude is related by crossing to meson photo- 
production £>7 — > ir°p. The latter process has been recently investigated within the handbag 
approach in Refs. [8, 9]. It turned out that a leading-twist calculation of the partonic sub- 
process, which is 7T° photoproduction off quarks, is insufficent. In leading- twist accuracy, 
i.e. considering only the one-gluon exchange mechanism for the generation of the meson, 
the resulting ^p — > 7T°p cross section is far too small. This parallels observations made in 
leading- twist calculations of the pion form factor [10]. The results are typically a factor 3 
to 4 below the admittedly poor data available at present. In fact the lowest-order Feyn- 
man graphs contributing to jp — > n°p within the handbag approach are the same as those 
occuring in the calculation of the pion form factor. Therefore, one has to conclude that for 
the existing data one is still far from the asymptotic region in which one-gluon exchange 
dominates and it is necessary to use a more general mechanism for the generation of the 
meson. Many alternatives and/or corrections to the leading-twist meson generation have 
been discussed in the literature reaching from higher-twist or power corrections to resum- 
mation of perturbative corrections. Still the description of meson generation is a matter 
of controversy. In order to remedy the situation a treatment of the subprocess is called 
for that does not postulate the dominance of any specific meson generation mechanism. 
Huang et al. [9] have proposed such a method. They utilized the covariant decomposition 
of the 79 — > n°q amplitudes proposed by Chew, Goldberger, Low and Nambu (CGLN) 
[11]. This decomposition separates the kinematic and helicity dependences of the subpro- 
cess from the dynamics of the meson generation which is embodied in the CGLN invariant 
functions. Exploiting properties of the invariant functions various dynamical mechanisms 
can be identified and consequences discussed. Comparison with experimental data on the 
ratio of n + and 7r~ photoproduction cross sections [12] provides evidence for the dominance 
of one invariant function out of the set of four. 

In this work we are going to investigate the corresponding time- like process pp — > 7 tt°. As 
for its space-like partner we will again make use of the CGLN decomposition as at least 
for the fixed-target option of FAIR (s < 30 GeV 2 ) but possibly also for the collider option 
(s < 210 GeV 2 ) the one-gluon exchange mechanism for meson generation will not dominate. 
Comparison with existing data from Fermilab [13] at 8.5 GeV 2 < s < 13.6 GeV 2 allows for 
a critical examination of the handbag charateristics. We will show that, in parallel to the 
space-like region, there are indications for the dominance of one of the invariant functions. 
In fact, this function is the s <-> t crossed one of the seemingly leading one in the space-like 
region. Morever, this function is the one that is fed by the leading-twist mechanism al- 
though not sufficiently strongly. Exploiting the crossing properties of the CGLN invariant 
functions, we are in the position to relate the space- and time-like processes quantitatively. 
Let us note that for the complementary kinematic regime, namely production of a virtual 
photon and a pion into the forward direction, a different factorization scheme was proposed 
by Pire and Szymanowski [14]. 
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2 The handbag amplitude for pp — > 7 7r 







The treatment of proton-antiproton annihilation into photon and meson parallels that of 
annihilation into two photons respectively its time-reversed process, two-photon annihila- 
tion into pp, investigated in [3]. Therefore, we can take over many of the results derived 
in that publication and for clarity we will also use a notation as close as possible to that 
employed in [3]. 

Obviously our crucial starting point is the assumption of handbag factorization of the am- 
plitude for the kinematical region s, —t, —u ^> A 2 where A is a typical hadronic scale of the 
order of 1 GeV. In this factorization scheme for which validity arguments have been given 
in Ref. [3], the process amplitudes appear as a product of a hard subprocess, q a q a — > 7 7r° 
and a soft pp — > qq transition matrix element which is parametrized by pp distribution 
amplitudes % = A,P,V,S introduced in [3] and where a is the flavour of the quark 
anti-quark pair emitted from the pp pair. The distribution amplitudes are time-like 
versions of general parton distributions for the proton. 

We are interested in large-angle scattering processes. For very high energies the handbag 
diagram will not dominate, but rather processes like that in Fig. la [15]. However, this 
is not a consequence of power counting, but rather of Sudakov suppression, which can 
only be expected to be effective at really large energies. (A nearly real parton entering or 
leaving a hard scattering process emits gluon-bremstrahlung. Excluding this part of the 
cross section by requireing exclusivity therefore leads to suppression.) See [16] for details. 
The asymptotic contribution a la Brodsky and Lepage [15] to our process has not yet been 
worked out. Experience with other processes like 77 — > pp [17] or ^p — > ^p [18] let us, 
however, expect that this contribution to the cross section is likely to lie way below the 
experimental points for s of the order of 10 GeV 2 . 

As we will follow closely the discussion of two-photon annihilation into pp, presented in 
[3], we refrain from presenting all details of the theoretical analysis, and rather restricting 
ourselves to a sketch of the main points. We work in a symmetric frame where the momenta 
of the proton and the antiproton are defined in light-cone coordinates as 



where f3 — y 1 — 4m 2 /s and m is the mass of the proton. In the following we will work in 
the massless limit. The transverse direction is characterized by the two-dimensional vector 
e_|_ for which we choose (1, 0). The important feature of a symmetric frame in the time-like 
region is that the skewness is 




(1) 





(2) 



The momenta of the photon and the pion read 
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Figure 1: a) One of the leading-twist contributions at asymptoticaly large s. The two blops 
on the left-hand side represent ordinary proton distribution amplitudes while the one on 
the right-hand side is the pion distribution amplitude, b) The handbag contribution at 
large but non-asymptotic s. The blob at the left-hand side represents the pp distribution 
amplitude, the one on the right-hand side the general dynamics for qq — > 771" parametrized 
in terms of the CGLN invariant functions. The momenta and helicities of the various 
particles are indicated. 



in our frame of reference which is also a cm. frame. The scattering angle is denoted by 6 
and the trigonometric functions of it are related to the Mandelstam variables by 



• n %VtU t-U 

sin U = , cos v = . (4) 

s s 

The starting point for the derivation of the handbag amplitude is the expression 

M = Y,ee a Jd*kJ ^e"^ ( | T*°(0) | p(p) p(p') > , (5) 

where we omitted helicity labels here for convenience. The hard scattering kernel, de- 
scribing the subprocess q a (k) q a (k') — > 771" , is denoted by 7i a . The sum runs over the 
quark flavours a = u,d, s; e a is the corresponding charge in units of the positron charge 
e. As discussed in detail in [2, 3, 19], the pp — > qq transition at large s can only be soft 
if the outgoing quark and antiquark have small virtualities and each carries approxima- 
tively the momentum of the baryon or antibaryon; the deviations of the parton momenta 
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from the hadronic ones in size and direction are of order A 2 /s. It can be further shown 
that the dominant Dirac structure of the soft transition matrix element in (5) involves 
the good components of the quark fields in the parlance of light-cone quantization. The 
hard scattering amplitudes can approximately be calculated with on-shell parton momenta: 
H a = u a H a v a . This guarantees gauge invariance. Within our calculational scheme we also 
use the approximation that the proton (antiproton) dominantly emits a valence quark (an- 
tiquark). Contributions from the emission of fast sea quarks are expected to be small. 
Putting all this together we obtain, in full analogy to the case of pp — > 77 [3], the handbag 
amplitude for the process under consideration (for the sake of legibility explicit helicities 
are labelled by their signs) 

AW = Y,l e «{{ H Zo, + - + H »o,- + ] S ^' F v*(s) (6) 

a z 

+ [#,V- - tf,V + ] [2^_,, (Fr(s) + F«*(s)) - £ 5 VV , Fr(s)} ) + O(^) , 



where we make use of the definition of the annihilation form factors in terms of the pp 
distribution amplitudes $f as introduced in [3] 

F t a (s) = f'dz^izX.s,^) for t = V,A,P. (7) 
Jo 

In Ref. [3] the form factors are defined for the soft transitions q a q a — > pp while here we 
consider the time- reversed transitions. Therefore, the complex conjugated form factors 
occur here. The fourth form factor or distribution amplitude, the scalar one, decouples in 
the symmetric frame (C = §)> see Pi- Note that the relations (7) hold for any physical 
value of the skewness. They also hold for any value of the factorization scale, /j,p, of the 
distribution amplitudes, since the vector and axial vector currents have zero anomalous 
dimensions. Because of symmetry properties of the pp distribution amplitudes [3], F£ and 
Fp project onto the C-even part of the pp state while Fy is related to the C-odd part. 
On the other hand, the 771" state is C-odd. This apparent conflict is a consequence of 
neglecting configurations where the antiquark is emitted in the direction of the proton 
instead of the antiproton. According to the arguments presented in [3] this contribution is 
small. 

We emphasize that quark-antiquark helicity non-flip does not occur in (6) although there 
is baryon helicity non-flip. Although we deal with massless quarks, quark helicity non-flip 
could contribute since we have not specified the mechanism that controls the generation 
of the meson. While the leading- twist, one-gluon exchange mechanism is pure helicity 
flip, twist-3 effects, for instance, produce non-flip contributions. Such contributions which 
are accompanied by helicity-non-flip generalized distribution amplitudes and associated 
form factors [20], have been studied for the corresponding space-like process in [9]. In 
this respect pp annihilation into 7 n° differs from annihilation into two photons [3] . The 
subprocess for the latter reaction is pure helicity flip for massless quarks. In parallel to pion 
photoproduction where quark helicity flip is neglected [9], we assume that quark helicity 
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non-flip contributions are negligible in the time-like region. In principle, this conjecture can 
be tested by measuring helicity correlations. As we argued above the dominant contribution 
to our process comes from the emission of valence (anti)quarks by the (anti)proton. We 
therefore have to take into account only the two subprocesses uu — > 7 tt° and dd — > 7 tt°. By 
isospin symmetry the weights of these two subprocesses are C" = I/a/2 and C^ = — 1/V2, 
see the discussion in [9]. It is therefore convenient to pull out these weight factors from 
the subprocess amplitudes, H a = C® H, and absorb them as well as the corresponding 
fractional charges into the form factors by introducing, following previous work [3, 8, 9, 19], 
annihilation form factors specific to pp annihilation into a photon and a 7r° (i = V, A, P) 



R 



r° 1 ( 



e F" 



(8) 



After this procedure the subprocess amplitudes are flavor independent and we can drop 
the corresponding superscript. 

To proceed we make use of the CGLN covariant decomposition of the subprocess am- 
plitudes [11]: 



H. 



J2C i (s,t)u(k : \)Q i v(k',\'), 



(9) 



i=l 



with the manifestly (elm.) gauge invariant covariants 

Qi = -275 k ■ e* q ■ q' -k ■ qq' ■ e* 
Q 2 = -2 75 \k.qf-k-e*i\ , 



Q 3 = 
Qa = 



75 



(10) 



encoding the helicity structure of the subprocess. The vector e denotes the polarisation of 
the photon and 



k = -(k — k') 



;n) 



The invariant functions Cj depend on the detailed, unknown dynamics. Using the spinor 
definitions given in Ref. [20], Eq. (17), one can calculate the helicity amplitudes in terms 
of these invariant functions. This reveals that C\ and C 4 only contribute to the parton 
non-flip amplitudes in the massless limit which, as just discussed, will be neglected. The 
other two invariant functions contribute to the helicity flip amplitudes 



H. 



+o,- 



H. 



+0, 



(s,t) 



(s,t) 



s 



Co — C? 



C2 + Cz 



;i2) 



The Mandelstam variables in the subprocess may differ from the ones for the full process 
due to the proton mass, see for instance [21]. Here, in this work we will ignore such possible 
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target mass corrections which are of order m 2 /s. Inserting (8) and (12) into (6), we obtain 
the amplitudes 



M 



{[+(t-u)C 2 -sC 3 



8 V - V l Ry 



+ 



sC 2 -(t-u)C 3 ]\2u <W (Rf * + Rf *) - ^ 5„ v , R? 

J 1 2m 



, (13) 



which lead to the pp — > 7 7T° cross section 



da 
~dt 



(pp -> 7 7T°) 



^clm 
~32~ 



X 



{\C,(s,t)\ 2 [(^) 2 \Rf(s)\ 2 + (R::(s)f 



+ |C 3 (M)| S 



\R V (s) 



t — U\2 



S 

- 2Re[C 2 (s,t)Cl(s,t)](^) [^^(^P + TO ^)) 2 ]} • (14) 

neglecting 0(m 2 / s) terms. Since always a particular combination of the form factors R^° 
and Rp appears in the cross section we introduce an effective form factor 



R cS - (\ R A + R P I 2 + 



4m 2 



|fl?T) 



1/2 



(15) 



Lack of suitable polarization data prevents to disentangle both these form factors anyway. 
The eqs. (13), (14) are the s <-> t crossed versions of the handbag amplitudes and the cross 
section for ^p — * 7r°p derived in Ref. [9]. There are only some minor modifications for 
the form factors occuring as a consequence of the specific frame of reference used in the 
space- and time-like regions. These modifications which have extensively been discussed 
in Ref. [3], are: Instead of vector and tensor generalized distribution amplitudes or form 
factors it is of advantage to apply the Gordon decomposition to the vector piece in the 
time-like region and trade the tensor form factor for the scalar one. It then turns out 
that, in contrast to the space- like region, the scalar form factor decouples in the symmetric 
frame with ( = 1/2 while the pseudoscalar one contributes. 

As is obvious from Eq. (6) the pseudoscalar anihilation form factor R P ° generates the 
pp — > qq transitions where proton and antiproton have the same helicities while quark 
and antiquark have opposite ones. This implies parton configurations of the pp system 
with non-zero orbital angular momentum. It is therefore expected that, at large s, the 
pseudoscalar form factor is suppressed as compared to the other form factors by at least 
l/y/s. If so it will not dominate over the other terms in Eqs. (14), (15) with increasing s. 

3 The annihilation form factors 



As is obvious from Eq. (13) an analyis of the process pp — > 7 tt° at large s, —t, —u requires 
information on the annihilation form factors (8) which encode the physics of the soft 
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pp — > qq transition. Presently the form factors cannot be calculated from first principles 
in QCD. Also there is no model calculation available for them currently. Thus, we have to 
determine the form factors phenomenologically. 

Exploiting the universality property of the generalized distribution amplitudes, we can use 
the form factors that have been determined from an analysis of two-photon annihilations 
into pp pairs within the handbag approach [3]. We however do not use the results of the 
numerical analysis presented in Ref. [3] since we have now at our disposal the very accurate 
BELLE data on the differential and integrated cross sections for 77 — > pp [7]. These data 
allow a determination of the vector and effective form factors from a Rosenbluth type 
separation. We therefore redo the analysis of the annihilation form factors. In the handbag 
approach [3] the 77 — > pp cross section is given by an expression which is analogous to the 
I C2 1 -term in Eq. (14) with form factors defined analogously to Eqs. (8) and (15): 

| (77 ^ pP) .l^_L_ {lW + Wco ^ } . (16) 

A fit to the BELLE data on the differential cross section at the highest measured energy 
(3 < y/s < 4 GeV , implying an average s of s = 10.4 GeV 2 ) and on the integrated cross 
sections for s > 8 GeV 2 provides 1 > 2 

s 2 R? s = (2.90 ±0.31) GeV 4 (s/s ) { - 1AOMA5) , 

\s 2 R^\ = (8.20 ±0.77) GeV 4 (s/s ) (_1 - 10±ai5) , (17) 

where the same energy dependence is assumed for both form factors. The fit to the dif- 
ferential cross section is compared to experiment in Fig. 2. The quality of the fit is, with 
the exception of the data point at cos6> = 0.55, very good. The failure with this point 
may be regarded as an indication that the handbag approach is here at its limits. Indeed 
for s ~ 10.4 GeV 2 this scattering angle corresponds to t ~ — 2 GeV 2 which is not much 
larger than A 2 . For comparison we also show in Fig. 2 the scattering angle dependence 
of the cross section in the case that R\ equals R v (this curve is arbitrarily normalized). 
Evidently the form factors are different. 

The BELLE data exhibit clear evidence for violations of the dimensional counting rules 
[22], the integrated cross section falls off faster than s~ 5 resulting in an energy dependence 
of the annihilation form factors as given in (17) instead of a constant behaviour. As we 
mentioned above, due to Sudakov suppressions the annihilation form factors are expected 
to decrease faster than 1/s 2 for very large s. Recent BELLE measurements [23] on 77 — > 
7r + 7r~ and K + K~ as well as the JLab measurement [24] of Compton scattering also signal 
violations of the dimensional counting rules for hard scales of the order of 10 GeV 2 . 
We note that the values for the form factors quoted in (17) are somewhat different from 
the estimate given in Ref. [3] . This estimate was based on the the CLEO and VENUS data 

1 The cross section data in the region of r\ c formation is removed from the fit. Possible signals from the 
P-wave charmonia are not visible in the BELLE data. 

2 In this energy region the BELLE data for the differential cross section reveal a minimum at an scat- 
tering angle of 90°. For lower energies the cross sections behave differently [7]. 



8 




Figure 2: The differential cross section for 77 — > pp versus cos# at 3.0 < y^i < 4.0 GeV. 
Data taken from Ref. [7]. The solid line represents our fit while the dashed one shows the 
angle dependence assuming form factors, R^ s = Ry, for comparison. 

[6] which are of markedly worse quality than the BELLE data [7] and were only available 
at the rather low value s = 7.3 GeV 2 . 

Since we have only protons at our disposal an exact flavor decomposition is not possible. 
We therefore follow Ref. [3] and simplify the expressions by taking a single real constant p 
for the d/u ratio of all proton form factors 

Ff = pF?, i = V,A,P. (18) 

This ansatz parallels the behavior of fragmentation functions for d — > p and u — > p transi- 
tions. In Ref. [3] the range of values 

p = 0.25 - 0.75 (19) 

has been considered. We remark that simple quark counting arguments give p = 1/2 [25]. 
Combining (18) with the flavor decomposition of the form factors for 77 — > pp 

R? = elFf + e'F*, (20) 

and neglecting non-valence quark contributions, one finds for the form factors (8) of 771" 
production 

Rf - t ] ~ ( ed iX i = V,A,P (21) 

e u V2 1 + {e d /e u ) 2 p 

Using the numerical values for the RJ quoted in (17) and adding the errors of the form 
factors and p quadratically, we arrive at the estimate 

Is 2 ^ ! = (3.42 ±0.4) GeV 4 (s/s )~ L10±a15 , 

|s 2 i^°| = (9.67 ±1.1) GeV 4 (s/s )" 110±015 . (22) 
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As we said above the effective form factor is likely dominated by the axial-vector form 
factor since the pseudoscalar one involves parton orbital angular momentum and is, there- 
fore, expected to be suppressed as compared with the other form factors at large s. An 
experimental separation of the axial vector and the pseudoscalar from factors requires po- 
larization experiments. The analysis of such observables is beyond the scope of this work. 
We will only briefly comment on this issue below. For the analysis of cross sections we 
only need the effective form factor. 

Additional information on the vector form factor is provided by the data on the magnetic 
proton form factor G V M in the time-like region for s up to 14.4 GeV 2 [26]. This form factor 
is related to the vector form factors (7) by 3 

G P M (s) = E ^F v a (s). (23) 

a=u,d 

The E835 data on the scaled magnetic form factor s 2 \Gm\ can also be represented by a 
power law 

|s 2 G^| = (2.46 ±0.16) GeV 4 ( S /s ) ( - a75±a34) . (24) 

An estimate of Ry from G^ along the same lines as described above, leads to a smaller 
value than given in (22) [3]. Given the approximations made in the derivation of the 
handbag contribution [3], there is however no real contradiction. 



4 Analysis of the Fermilab (E760) data 

Let us turn now to the discussion of the CGLN invariant functions. They have definite 
behaviour under t <-> u crossing [11]: 

C 2 {s,t) = C 2 (s,u), C 3 (s,t) = -C 3 {s,u), (25) 

The invariant functions have dynamical singularities. This can be seen from their leading- 
twist contributions which are expected to dominate at large s. In this kinematical region 
the subprocess qq — > 771° is dominated by one-gluon exchange to be calculated in collinear 
approximation to lowest order of perturbative QCD. One obtains 

C c 2 °\s,t) = ^, C c 3 °\s,t) = 0, (26) 

where 

a 2 co11 = 4W/4)/^UAV (27) 

Here, n 2 R is an appropriate renormalization scale, f n the pion decay constant, Cp = 
(iV 2 - 1)/(2N C ) is the usual SU(N C = 3) color factor and the last factor in (27) is the 

3 Note that only the form factors (7) are universal in contrast to the flavor combinations occuring in 
specific processes. 
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1/r moment of the pion distribution amplitude. To leading-twist accuracy also C\ and C4 
are zero. These invariant functions control the situation where quark and antiquark have 
the same helicity. They are only non-zero at twist-3 (or higher) level. As we mentioned 
in the introduction the leading-twist, lowest order pQCD contribution to the subprocess 
falls short as compared to experiment. It is therefore suggestive to assume that handbag 
factorization holds and that a more general mechanism than leading- twist, unknown at 
present, is at work for the generation of the meson and enhances the invariant function C 2 
sufficiently 

In view of this let us try now to find a suitable parametrization of the invariant functions. 
Their singularity structure will not be altered by the inclusion of more complicated dy- 
namical effects such as higher-orders of perturbative QCD, transverse degrees of freedom 
and/or higher twists or power corrections. For instance, by insertion of an infinite number 
of fermionic loops in the hard gluon propagator and by interpreting the ambiguities in 
the resummation of these loop effects - known as infrared renormalons - as a model of 
higher-twist contributions an enhancement of the leading-twist by a large factor may be 
obtained [27]. Therefore, a general ansatz for C 2 reads 

C 2 (s,t) = ^f 2 (t,u), (28) 
tu 

where f 2 is a symmetric function of t and u. It may comprise powers of tu/s 2 , i.e. powers 
of sin#, or symmetric combinations of hit and lnw. 

From the leading-twist result (26) one may expect that the invariant function C 3 plays 
only a minor role at s of the order of 10 GeV 2 and may be neglected. This assertion is 
further supported by an observation made for 7r ± photoproduction [9]: the experimental 
ratio of these two processes is, within the handbag approach, only understood if \C 2 \ 
I C 3 1 provided quark helicity flip can be neglected. Here, the are the s <-> t crossed 
invariant functions Cj. With regard to this observation we will first discuss a scenario with 
a dominant C 2 assuming for it the ansatz (28) with the simplest choice f 2 = 1. Taking the 
estimate (22) for the annihilation form factors we are left with only one free parameter, 
namely the normalization constant a 2 . In general a 2 is a complex number but the cross 
section only probes its modulus. The energy and angle dependence of the pp — > 771- cross 
section is fixed by the handbag physics. 

The differential and integrated cross section read for this scenario (9 is the angle between 
incoming proton and the outgoing photon in the center of momentum frame.) 



da a e \ m \a 2 



2 



dcos6 ' 4s 6 sin 4 6» 

|2 



| S 2 <;| 2 + cos 2 ^|s 2 <| 2 



a = 



Ctclm \ a 2\ 

1-4 



iln^(| S 2 ^| 2 -| S 2 <| 2 ) 

(29) 



+ ^(\s 2 R^\ 2 + \s 2 R* v °\ 2 
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s[GeV] 



Figure 3: The integrated cross section for pp — > 771"° versus s. Filled (open) circles: 
integrated over the range | cos^l = to z , z = 0.5 (0.6). Data taken from [13]. The black 
line represents the prediction from the handbag approach with the error bands evaluated 
from the uncertainties of the annihilation from factors. The dashed line represents a power- 
law fit for z = 0.6 with a ~ (s/s ) 8 - 16±a12 . 



where zq = | cos(9o| is the limit of integration. With zq = 0.5 the integrated cross section is 

12 



a = 864nb GeV 2 



\a 2 \ 



R* s \ 2 + 0.096 \s 2 Rl 



(30) 



A fit to the E760 cross section data provides the value (13.39±0.10) GeV for the parameter 
1 0-2 1 - The fit is compared to the data [13] in Fig. 3. With regard to the errors of the 
annihilation from factors quoted in (22), very good agreement of theory and experiment 
can be claimed. The contribution from J/\& formation is very small and can be ignored. 
We also neglect a possible contribution from the Its decay into 771"° has not yet been 
observed. 

For a few energies one may also integrate the E760 data up to z = 0.6. For comparison 
and in order to demonstrate the internal consistency of our numerical analysis we fit a 
power law to these data. As can be seen from Fig. 3 the fit is in good agreement with 
the data and yields the value of 8.16 ± 0.12 for the power in perfect agreement with the 
handbag approach and the form factors (22) estimated from 77 — > pp. We note in passing 
that a power of 6 is expected from the dimensional counting rules our process. 
The results from the handbag approach for the differential cross sections are shown for two 
energies and for | cos#| < 0.6 in Fig. 4. The excellent agreement between the theoretical 
results and the E760 data is obvious. Similarly good results are obtained for all energies 
above 10 GeV 2 while for lower energies the cross section data exhibit marked fluctuations 
which might be indicative of prominent low orbital angular momentum partial waves. This 
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Figure 4: The differential cross section for pp — > 771" versus cos 9 at s = 12.44 GeV 2 
(left) and 13.59 GeV 2 (right). Data taken from [13]. The solid lines with the error bands 
represent the prediction from the handbag approach. For comparison results are also shown 
where the cross section behaves oc 1/ sin 2 ^ instead of oc 1/ sin 4 ^ (dashed line). 



behaviour of the pp — > 771" cross section parallels that observed in pp — > 77 [6, 7] and is the 
reason why we determined the vector form factor only from BELLE data [7] for the cross 
section at s = 10.4 GeV 2 . In order to elucidate the structure of the handbag contribution 
further we display the pp — > 771" cross section at a scattering angle of 90° in Fig. 5. Also 
here we observe fair agreement between the handbag approach and the data although the 
errors are substantial. 

In Fig. 4 we display for comparison also the results of a calculation assuming f 2 = sin 9 
( for simplicity we retain the value of the parameter a 2 ). In this case the angle dependence 
of the cross section coincides with that for pp — > 77 but it seems to be too weak for the 
77T channel in contrast to the two-photon channel where it is in very good agreement with 
the BELLE data for | cos#| < 0.5. Thus, one may conclude that, within errors, there is no 
need for f 2 7^ 1. In other words, an ansatz for the invariant function C 2 that resembles the 
singularity structure of the leading-twist result but with a strength adjusted to experiment, 
describes the data very well. 

From Fig. 5 it is also obvious that a scenario | C* 3 1 ^> \C 2 \ is in conflict with experiment. The 
zero of the t <-> u antisymmetric invariant function C3 at 90° would lead to a corresponding 
zero in the differential cross section which is not seen with the exception of a few energies 
where a zero is possible within errors 4 . It goes without saying that a small admixture 
of C3 to a dominant C 2 cannot be ruled out. There is yet another test of the internal 
consistency of the handbag approach. Our process is related to photoproduction of pions 

4 Explicit fits of the C3 scenario to the differential cross section data, using for instance the ansatz 
C3 = 53 cos 6 /(tu), provide results which are too small in the vicinity of 90° as a consequence of this zero. 
On the other hand, the fit to the integrated cross section is of the same quality as in the C 2 scenario (with 
|o 3 | = 20.05 ±0.15 GeV). 
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Figure 5: The differential cross section for pp — ■> 7 tt° at 90° versus s. Data, representing 
the average of the cross sections in the two bins adjacent to 90°, are taken from [13]. The 
solid line with the error band represents the prediction from the handbag approach. 

by s <-> t crossing as we mentioned occasionally. The form factors are functions of t in 
the latter case and are known from a recent analysis of nucleon form factors exploiting the 
sum rules satisfied by generalized parton distributions [16]. Using the invariant function 
related by s <-> t crossing to C2 

C a (M) = — , (31) 
su 

where the Mandelstam variables are now those for photoproduction, one finds fair agree- 
ment with the high-energy wide-angle SLAC data [28] if the parameter a 2 is adjusted to 
these data (a 2 = 20.3 GeV). For very large s the parameter a 2 should coincide with 
1 0,2 1 = 13.39 ±0.10 GeV but at the actual s values of order 10 GeV 2 this discrepancy is not 
implausible. Typical differences between time- and space-like values of many quantities 
are of this size [29]. It is, however, important to realize that the photoprodcution data [28] 
need confirmation. 



5 Predictions for FAIR/GSI 

Obviously Figs. 3 and 5 provide already predictions for the cross sections to be expected at 
FAIR. In Fig. 6 our prediction for the differential cross section at s = 20 GeV 2 is shown (the 
fixed angle cross section scales as s 8 - 2±0 - 3 in our approach). As the discussed luminosities 
reach up to 2- 10 4 nb _1 /day for unpolarized reactions and 10 3 nb _1 /day for polarized ones, 
the data situation will improve drastically, ones FAIR is operational [5]. 
The fact that the protons and antiprotons will be polarized adds very attractive additional 
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Figure 6: Cross section prediction for the FAIR fixed-target energy s = 20 GeV 2 . 

possibilities. It allows e.g to separate and R P ° . In anlogy to the two-photon channel 
[3] the helicity correlation All between proton and antiproton is given by 

da(++) - da{+-) 
LL ~ da(++) + da(+-) 
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(32) 



where da{vv') is the cross section for polarized proton-antiproton annihilation. 
Our analysis can be easily generalized to other pseudoscalar mesons, e.g. to the 7 r\ (7/) 
channel. The only complication arises from the additional ss and two-gluon Fock com- 
ponent the 77 and 77' mesons possess. For partons which are emitted nearly collinearly by 
the proton and antiproton, each carrying a large fraction of its parent's momentum, the 
strange quarks are likely strongly suppressed. Not much is known about the two-gluon 
component. There is only some, not very precise information on its leading-twist distri- 
bution amplitude from a next-to-leading order analysis of the 7 — > 77, 77' transition form 
factors and the inclusive T — > r)'X decay [30]. These analyses tell us that to leading- twist 
accuracy the contribution from the two-gluon Fock component to v( production may be 
sizeable while it is strongly suppressed in the case of the 77. Beyond the leading-twist level 
the role of the two-gluon contribution is unknown in the case of the 77'. Since the 77 meson 
is dominantly a flavor octet state, however, its two-gluon component should be suppressed 
in any case. In view of this one can treat the rj meson in the quark flavor basis and, 
neglecting the strange component, regard it as cos0 r] q where is the usual mixing angle 
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in the quark-flavor basis and r] q is the quark part of the 77 wave function. The annihilation 
form factors for the production of the isoscalar r\ q read 

R? = ^(e u F t u +e d F t d ). (33) 

Their numerical values are somewhat smaller than those for the 7r° given in (22). Up to 
this difference and a new value of the parameter a 2 the handbag approach then predicts 
the same energy and angle dependence of the cross sections for the 777 channel as for the 
77T° one. Provided the two-gluon component also plays only a minor role in the case of the 
rj the ratio of the 777 and 777' cross sections is given by 

dajpp -> 777) = c 
da(pp — » 777') 

In [31] the 77 — 77' mixing angle has been determined to be 39.3°. 

Our results for the •yn channel can also be straightforwardedly generalized to the 7V1, 
channel where Vl is a longitudinally polarized vector meson, if we rely again on valence 
quark dominance. The annihilation form factors for p° and uj production are the same 
as given in (22) and (33), respectively. In contrast <fi meson production is expected to be 
strongly suppressed because of the mismatch of the proton and meson valence quarks. 
The cross section for the ^Vl channel reads 



da VL Ctelm \d2 L \ 2 



dcosO 4s 6 sin 4 ^ 



|s 2 i^| 2 + cos 2 0| S 2 <|| 2 l. (35) 



The dependencies on the vector and the effective form factors are reversed in this case 
as a consequence of parity invariance (see, for instance Ref. [8]). We therefore expect a 
somewhat flatter angular dependence, close to 1/ sin 4 9, for the vector meson channels than 
for 77T°. The isolation of longitudinally polarized p° mesons might be possible with the 
planned PANDA detector [5]. 

The generalization to transversally polarized vector mesons, Vt, is more intricate. One 
has to consider the subprocesses qq — > yVr with equal and opposite quark and antiquark 
helicities. In the first case one has to introduce new pp distribution amplitudes and, hence, 
a new set of associated annihilation from factors. These new distribution amplitudes are 
time-like versions of the helicity-flip GPDs introduced in Ref. [20]. For opposite quark and 
antiquark helicities, on the other hand, the formation of the transversally polarized vector 
meson requires a different mechanism than for longitudinally ones. It is beyond the scope 
of this work to analyse the processes pp — > yVr- 



6 Summary 

We have analysed the reaction pp — > 771"° for large scattering angles assuming handbag- 
dominance. We obtained a rather satisfying description of the existing Fermilab data from 
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E760. Far more precise data can be expected from FAIR/GSI which also provides the 
possibility to use polarization variables to separate the different contributions. Such exper- 
iments would contribute to the determination of the generalized distribution amplitudes of 
the nucleon. It would also allow to clarify the nature of the dominant reaction mechanisms 
as a function of cm energy s. We argue that a similar analysis can be performed without 
major problems for annihilation in a photon plus either another pseudoscalar meson or a 
longitudinal vector meson. 
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